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ABSTRACT 

We present results from a mini-survey of relatively high redshift (1.7 < z < 4) gravitationally lensed 
radio-quiet quasars observed with the Chandra X-ray Observatory and with XMM-Newton. The lensing 
magnification effect allows us to search for changes in quasar spectroscopic and flux variability properties 
with redshift over three orders of magnitude in intrinsic X-ray luminosity. It extends the study of 
quasar properties to unlensed X-ray flux levels as low as a few times 10~^^erg cm~^ s~^ in the observed 
0.4-8 keV band. For the first time, these observations of lensed quasars have provided medium to 
high signal-to-noise ratio X-ray spectra of a sample of relatively high-redshift and low X-ray luminosity 
quasars. We find a possible correlation between the X-ray powerlaw photon index and X-ray luminosity 
of the gravitationally lensed radio-quiet quasar sample. The X-ray spectral slope steepens as the X-ray 
luminosity increases. This correlation is still significant when we combine our data with other samples 
of radio-quiet quasars with z > 1.5, especially in the low luminosity range between lO^'^-lO'*^'^ erg s~^. 
This result is surprising considering that such a correlation is not found for quasars with redshifts below 
1.5. We suggest that this correlation can be understood in the context of the hot-corona model for 
X-ray emission from quasar accretion disks, under the hypothesis that the quasars in our sample accrete 
very close to their Eddington limits and the observed luminosity range is set by the range of black 
hole masses (this hypothesis is consistent with recent predictions of semi-analytic models for quasar 
evolution). The upper limits of X-ray variability of our relatively high redshift sample of lensed quasars 
are consistent with the known correlation between variability and luminosity observed in Seyfert Is when 
this correlation is extrapolated to the larger luminosities of our sample. 



1. INTRODUCTION 

It is important to extend the study of quasars to high 
redshifts in order to understand the evolution of quasars 
and their environments. One of the main results of recent 
studies of quasar evolution is that the quasar luminosity 
function evolves strongly with redshift. (e.g. Boyle et al. 
1987). Many studies of quasar evolution are aimed at ex- 
plaining this luminosity evolution. The X-ray band probes 
the innermost regions of the central engine of the Active 
Galactic Nuclei (AGN). The study of AGN in X-rays may 
possibly answer the question of whether there is an evolu- 
tion in their central engines and how this is related to the 
evolution of the quasar luminosity function. The observed 
X-ray continuum emission of AGN is generally modeled 
by a power law of the form N{E) — Nq{E / Eq)"^ where 
N{E) is the number of photons per unit energy interval. 
Extensive studies of AGN during the past decade indicate 
that this power-law component is produced by Compton 
scattering of soft photons by hot electrons in a corona (e.g., 
Haardt & Maraschi 1993; Haardt, Maraschi, & Ghisellini 
1994). The study of this power-law component, its corre- 
lations with other AGN parameters, and its evolution re- 
veals important information on the accretion process of the 
central object. The X-ray photon indices of low-redshift 
radio-quiet quasars are measured to have mean values of 
~2.6-2.7 in the ROSAT soft X-ray band (Laor et al. 1997; 
Yuan et al. 1998) and -1.9-2.0 for the ASCA hard X-ray 
band (George et al. 2000; Reeves & Turner 2000), and are 
correlated with the H/3 FWHM (Reeves & Turner 2000). 
There is no strong evidence that the X-ray power-law index 
evolves with redshift or correlates with X-ray luminosity to 



date (George et al. 2000; Reeves & Turner 2000). Another 
important parameter that describes the broad band spec- 
tral shape of quasars is the optical-to-X-ray spectral index, 

quantified as ofox = log(/2kev//25Qo^)/ log(i'2kov/i^25oo 1^' 
where fauov and f ° are the flux densities at 2 keV 

jzKcv ■'2500 A 

and 2500 A in the quasar rest-frame, respectively. Re- 
cently, several studies have provided estimates of aox for 
very high redshift (z > 4) quasars (Vignali, Brandt, & 
Schneider 2003; Vignah et al. 2003a,b; Bechtold et al. 
2003). Vignah, Brandt, & Schneider (2003); Vignah et al. 
(2003a, b) found that aox is mainly dependent on the ultra- 
violet luminosity while Bechtold et al. (2003) found that 
Q^ox primarily evolves with redshift. In addition to spectral 
studies, variability studies of Seyfert galaxies show that the 
variability amplitude (excess variance) is anti-correlated 
with X-ray luminosity (Nandra et al. 1997; Leighly 1999). 
Variability studies have been extended to quasars by sev- 
eral groups (George et al. 2000; Almaini et al. 2000; Man- 
ners, Almaini, & Lawrence 2002). Low redshift quasars 
[z < 2) are found to have an excess variance consistent 
with the luminosity relation found in Seyfert Is and there 
is a possible upturn of X-ray variability for high redshift 
quasars with z > 2 (Almaini et al. 2000; Manners, Al- 
maini, & Lawrence 2002). 

It is also important to compare the properties of quasars 
near the peak of their comoving number density, thought 
to have occurred at z ~ 2, with low redshift quasars. This 
comparison may provide clues as to what caused the dra- 
matic decay of the quasar number density as the Universe 
expanded. 

Most of the observational and analysis techniques em- 
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ployed to date to study the evolution and emission mech- 
anism of faint, high-redshift quasars are based on ei- 
ther summing the individual spectra of many faint X-ray 
sources taken from a large and complete sample or obtain- 
ing deep X-ray observations of a few quasars. Although 
these techniques may yield important constraints on the 
average properties of high redshift quasars they each have 
significant limitations. 

Gravitational lensing provides an additional method for 
studying higli-rcdshift quasars. The extra flux magnifica- 
tion, from a few to ~ 100, provided by the lensing effect en- 
ables us to obtain high signal-to- noise ratio (S/N) spectra 
and light-curves of distant quasars with less observing time 
and allows us to search for changes in quasar spectroscopic 
properties and X-ray flux variability over three orders of 
magnitude in intrinsic X-ray luminosity. With the aid of 
lensing, we can probe lower flux levels than other flux lim- 
ited samples with similar instruments and exposures. This 
could be an important factor because the X-ray properties 
of high-redshift, low- luminosity quasars could be differ- 
ent from those of other quasars and by studying them we 
could possibly obtain information about the evolution of 
quasars. 

Similar lensing studies have also been performed in 
the sub-millimeter and CO bands (Barvainis & Ivison 
2002; Barvainis, Alloin, & Bremer 2002), which suggested 
that gravitational lensing could be an efficient method 
for studying various properties of high-redshift quasars. 
Here we present the results of an X-ray mini-survey of 
relatively high redshift gravitationally lensed radio-quiet 
quasars. We chose only the radio-quiet quasars in our 
sample because the powerful relativistic jets in the radio- 
loud quasars will introduce additional complication when 
modeling the continuum X-ray emission from the accretion 
disc. 

We use a Ho = 50 km s~^ Mpc~^ and go = 0.5 cosmol- 
ogy throughout the paper. 

2. OBSERVATIONS AND DATA REDUCTION 

Our mini-survey contains eleven gravitationally-lensed, 
radio-quiet quasars with redshifts ranging between 1.695 
and 3.911. Five of them contain Broad Absorption Lines 
(BALs) or mini-BALs in their rest-frame ultraviolet spec- 
tra. Most of the lensed quasars of our sample were ob- 
served with the Advanced CCD Imaging Spectrometer 
(ACIS; Garmire et al. 2003) onboard Chandra as part of 
a Guaranteed Time Observing program (Principal Inves- 
tigator: G. Garmire). The data for two of the them were 
obtained through the public Chandra archive. Several of 
them were observed twice. Three of the lensed quasars 
were also observed with XMM-Newton. Table 1 presents 
a log of observations, including redshifts, Galactic column 
densities, and exposure times. Each observation was per- 
formed continuously with no interruptions. 

All of the sources observed with Chandra were placed 
near the aim point of the ACIS-S array, which is on the 
back-illuminated S3 chip. All of the data were taken in 
FAINT or VERY FAINT mode. The Chandra data were 
reduced with the CI AO 2.3 software tools provided by 
the Chandra X-Ray Center (CXC) following the standard 
threads on the CXC website. Only photons with standard 
ASCA grades of 0, 2, 3, 4, 6 were used in the analysis. 



We used events in the 0.4-8 keV energy range in the spec- 
tral analysis and events in the 0.2-10 keV range for the 
variability studies. The source events were extracted from 
circles with radii ranging from 3" to 5" depending on the 
separations of the lensed images of each quasar. The cir- 
cles included all of the lensed images. Background events 
were extracted from annuli with inner and outer radii of 
10" and 30", respectively, centered on the sources. We 
adjusted the inner and outer radii of the background sub- 
traction annuli in some cases to avoid other sources in the 
fleld. The background contributes an insignificant amount 
to the count rate in a source region, even during a back- 
ground flare. The average exposure time for the Chandra 
observations was about 25 ks. The detected source count 
rates ranged between 0.003-0.08 s~^. 

Three targets were observed with the European Photon 
Imaging Camera PN and MOS detectors (Striider et al. 
2001; Turner et al. 2001) onboard XMM-NeAvton. The 
XMM-Newton data were analyzed with the standard anal- 
ysis software, SAS 5.3. The tasks epchain and emchain 
from SAS were used to reduce the PN and MOS data and 
photons of patterns < 4 and < 12 were selected from the 
PN and MOS data, respectively. The XMM-Newton data 
are affected more than the Chandra data by background 
flares because the Point Spread Function of XMM-Newton 
is significantly larger than that of Chandra- ACIS. Sev- 
eral strong background flares occurred during the XMM- 
Newton observations. These flares are filtered out in the 
spetral analysis. 

3. SPECTRAL ANALYSIS 

3.1. Power-law Continuum 

We performed spectral fitting in order to obtain the 
power-law indices of the X-ray continuum emission com- 
ponents of the quasars in our sample. The photon 
indices that we measured were in the observed 0.4-8 
keV energy range. We used the Chandra data for all 
of the lensed quasars in the spectral analysis except 
PG 1115-1-080, where we used the XMM-Newton data. 
RX J0911.4+0551 and APM 08279+5255 were also ob- 
served with XMM-Newton. We did not use these data, 
however, because there were large amplitude and long 
duration background fiares in the XMM-Newton obser- 
vations. The background flare in the XMM-Newton ob- 
servation of RX J0911. 4+0551 almost spans the entire 
observation. The background flare in the XMM-Newton 
observation of APM 08279+5255 occurs during the last 
^ 20 ks of the observation. We filtered the background 
flare time from this observation and performed a spec- 
tral analysis to compare with our Chandra results. For 
our later correlation analysis, we used the Chandra data 
of APM 08279+5255 to avoid the possible complication 
from the cross calibration between Chandra and XMM- 
Newton instruments. For those sources observed twice 
with Chandra, wc performed simultaneous fits to the spec- 
tra extracted from the two observations except in the case 
of APM 08279+5255, where we used only the second ob- 
servation because it was much longer than the first one. 

We extracted spectra for each quasar using the CIAO 
tool psextract. We extracted events from all of the lensed 
images for each source except for the case of H 1413+117. 
A possible microlensing event in H 1413+117 appears to 
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amplify and distort the spectrum of image A only and 
thus affects the spectral slope greatly (Chartas et al. 
2003). In this case, we extracted the spectrum of the mi- 
crolensed image A and the spectrum of the other three 
images separately. A microlensing event is also detected 
in Q 2237+0305. However, the spectral shape of the mi- 
crolensed image A of Q 2237+0305 is not significantly af- 
fected by the microlensing event (Dai et al. 2003). 

Spectral fitting was performed with XSPEC VI 1.2 (Ar- 
naud 1996). There arc typically several lumdreds (180- 
6000) detected source events in the spectra of the target 
quasars, except for the spectrum of HE 2149—2745 which 
has only 23 detected events. The moderate S/N of our 
spectra allows us to fit each of them individually using 
relatively complex models. Thus we can constrain the un- 
derlying power-law slopes more accurately than in previ- 
ous studies of unlensed quasars of similar redshift. The 
models we used arc listed in Column 3 of Tabic 2. All 
spectral models included an underlying power-law model 
modified by Galactic absorption. The Galactic absorbing 
columns were obtained from Dickey & Lockman (1990). 

To account for the recently observed quantum efficiency 
decay of ACIS, possibly caused by molecular contamina- 
tion of the ACIS filters, we have applied a time-dependent 
correction to the ACIS quantum eflaciency implemented 
in the XSPEC model ACISABSl.l.^ If the source was a 
BALQSO with a medium S/N spectrum, we added a neu- 
tral absorption component at the redshift of the source. 
For high S/N BALQSO spectra such as PG 1115+080 and 
APM 08279+5255 we used the absori model in XSPEC to 
model the intrinsic absorption component as an ionized 
absorber (Chartas ct al. 2002; Chartas, Brandt, & Gal- 
lagher 2003). Wc added a neutral absorption component 
at the redshift of the lens for Q 2237+0305 (Dai et al. 
2003). For some of the spectra containing emission or ab- 
sorption line features, we added Gaussian line components 
to model them accordingly. The spectra of the quasars 
are displayed in Figures 1. The spectral fitting results are 
given in Table 2. In Table 2, we also list the observed 0.4-8 
keV fluxes and rest-frame 0.2-2 keV and 2-10 keV lumi- 
nosities for the lensed quasars in the sample not corrected 
for the lensing magnification. The unlensed luminosities 
are calculated in §4. We have corrected for the various ab- 
sorption components when calculating X-ray luminosities. 
These absorption components include Galactic absorption, 
absorption in the k;nsing galaxy, intrinsic absorption, and 
the ACIS contamination. 

The mean photon index of the radio-quiet quasars in 
our sample is 1.78 ± 0.06 with a dispersion of 0.27 ± 0.08 
and the median photon index of the radio-quiet quasars in 
our sample is 1.86. 

3.2. aox 

We calculated the optical-to- X-ray power-law slope, aox, 
for the quasars in our sample. The differential magnifica- 
tion between the optical and X-ray band is insignificant 
because both of the source regions are estimated to be 
much smaller than the Einstein radius. The rest-frame 2 



keV flux densities were calculated from the best-flt models 
to the spectra of the quasars. The redshifts of our sam- 
ple range from 1.7 to 4, thus the rest-frame 2 keV energy 
falls in the 0.4-0.74 keV observed-frame energy range. We 
removed all of the absorption components including the 
intrinsic ones for the BALQSOs when calculating the rest- 
frame 2 keV flux densities in order to study the properties 
of the intrinsic continuum. The rest-frame 2500 A fluxes 
were obtained from the published optical data in the liter- 
ature and from current ongoing programs such as the CfA- 
Arizona Space Telescope LEns Survey (CASTLES) ^ and 
the Optical Gravitational Lensing Experiment (OGLE). ^ 
We first converted the optical magnitudes to flux densities 
at the effective wavelengths of the filters, then extrapo- 
lated them to the rest-frame 2500 A flux densities. We 
used the standard relations between magnitudes and flux 
densities for the V and R bands. For the F814W mag- 
nitudes obtained with the HST WFPC2, we used the re- 
lation from Holtzman et al. (1995) to convert between 
magnitudes and flux densities. We used oc z^"" "^ in 
the extrapolation. Typical quasar optical spectral indices 
are in the -0.5 to -0.9 range (Schneider et al. 2001). We 
chose the optical magnitudes from the bands closest to 
the redshifted 2500 A wavelength to reduce the error in 
the extrapolation. We corrected the Galactic extinction 
based on Schlegel, Finkbeiner, & Davis (1998). The ex- 
tinction induced by the dust from the lensing galaxies is 
not significant except for Q 2237+0305, which is lensed by 
the central bulge of the galaxy. We corrected it according 
to the appendix of Agol, Jones, & Blaes (2000). The aox 
values are listed in Table 3. The listed error-bars of oiox are 
at the 1(7 level and account for the uncertainties in the X- 
ray spectral slope and normalization, the conversion from 
optical magnitude to flux density and the extrapolation 
of the optical/UV spectrum. We note that a source may 
have varied between the optical and X-ray observations 
thus leading to systematic errors in our estimates of aox- 
For the range of aox in our sample, a change of optical 
flux by a factor of two will lead to a change of Aaox^ 0.1. 
The rest-frame 2 keV and 2500A flux densities and rest- 
frame 2500A luminosity densities are also listed in Table 3. 
These values are corrected for lensing magniflcation as we 
discuss further in §4. In conclusion we find that the mean 
aox value for all the quasars in our sample is —1.70 + 0.02 
(-1.66 ± 0.02 without HE 2149-2745 which has only 23 
X-ray events) with a dispersion of 0.20 ± 0.03. 

4. MAGNIFICATION AND INTRINSIC X-RAY LUMINOSITIES 

We estimated the magnification for each of the lensed 

systems in order to obtain the unlensed X-ray luminosi- 
ties of the quasars. We first searched in the literature for 
magnifications of well-modeled systems. However, not all 
of the systems in our sample have been previously mod- 
eled in detail and for some cases the magnification values 
are not included in the published analysis. Wc modeled 
these systems with the gravlens 1.04 software tool devel- 
oped by C. Keeton (Kceton 2001a, b). We used a singular 
isothermal elliptical (SIE) mass profile with external shear 



^ ACISABS is an XSPEC model contributed to the Chandra users software exchange web-site http://asc.harvard.edu/cgi-gen/cont-soft/soft- 
hst.cgi. 

The CASTLES website is located at http://cfa-www.harvard.edu/glensdata/. 
^ The OGLE website is located at http://bulge.princeton.edu/''ogle. 
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in most cases, and took into account optical image posi- 
tions and flux ratios obtained from the CASTLES survey. 
We assumed a 20% error-bar for the flux ratios to account 
for uncertainties from microlensing or differences between 
optical and X-ray flux ratios. Our modeling results arc 
listed in Table 4. We compared the magnification values 
that wc obtained with those published in Barvainis & Ivi- 
son (2002) and found them to be in good agreement. We 
note that the magnification values obtained may have large 
systematic errors because the SIE model adopted in this 
analysis may not be suitable for all the systems. For ex- 
ample, in PG 1115-1-080, a set of different mass potentials 
are used in the modeling of the magnification in Impey 
et al. (1998). A magnification range of 20-46 was ob- 
tained for this system. Based on this example, wc adopt 
a factor of two as a typical systematic uncertainty in the 
magnification. In §6.1 we investigate how this system- 
atic uncertainty could affect for our conclusions via Monte 
Carlo simulations. 

5. VARIABILITY ANALYSIS 

The light-curves for the Chandra observations of the 
lensed quasars of our sample are displayed in Figure 2 
binned with a bin size of 1000 s (observed- frame). We 
show the light-curves of the longer Chandra exposure if 
two Chandra observations of the same target arc avail- 
able. We excluded HE 2149—2745 in our timing analysis 
because of its low S/N (only 23 events). 

To estimate the relative variability of the light-curves, 
we used the normalized excess variance (Nandra et al. 
1997; Turner et al. 1999; Edelson et al. 2002) defined 
as 

1 ^ 

i=l 

where N is the number of bins in the light-curve, Xi are 
the count rates per bin with error cTj, and jj, is the mean 
count rate. The error on a^^^ is given by S'D/[/i^(A'')^/^], 
where 

1 ^ 

^^ = Ar^E{[(^^-/^)'-^i]-^r-M'}' (2) 

i=l 

This method normalizes the variability amplitude to the 
flux, thus is less biased towards the low flux light- curves. 
It is better to apply this method to a sample of light-curves 
of similar lengths to avoid possible bias. For example, it 
would be easier to detect long time-scale variability from 
a light-curve with a long duration. In addition, the bin 
sizes should be similar in order to compare variability on 
similar time-scales. High S/N light-curves are needed for 
this analysis as the variability could be easily dominated 
by noise in the low S/N light-curves. This excess variance 
method also has some limitations. One of them is that it 
is not very sensitive to short flares of moderate amplitude. 
The excess variance values for Chandra light-curves with 
three different bin sizes are listed in Table 5. Wc also cal- 
culated the excess variance for the three sources observed 
with XMM-Newton and the results are listed in Table 6. 
The background was subtracted in the calculation. Over- 
all, the error-bars on the excess variance are quite large 
and in some cases, the values are negative. Therefore, for 
most quasars we can really only obtain an upper limit on 
the value of the excess variance. 



6. RESULTS AND DISCUSSION 

6.1. Luminosity and Spectral Index 

The unlensed 0.2-2 keV and 2-10 keV X-ray luminosi- 
ties of the lensed quasars in our mini-survey range from 
lO"*"^ to 10^^ erg s""*^. The mean spectral index of our sam- 
ple, 1.78, is consistent with the value of 1.84 from the 
recent study of very high redshift quasars (Vignali et al. 
2003b) and the value of 1.89 from the radio-quiet ASCA 
sample (Reeves & Turner 2000), especially when one con- 
siders the large dispersion of spectral indices in these stud- 
ies and our small sample size. 

We plot the rest-frame 0.2-2 keV and 2-10 keV X-ray 
luminosity against redshift in Figure 3(a), and the photon 
indices vs. rc;dshift, 0.2 -2 keV luminosity, and 2 10 keV 
luminosity in Figures 3(b), (c), and (d), respectively. Fig- 
ures 3(c) and (d) show a correlation between the spectral 
indices of our sample of lensed quasars and their X-ray 
luminosities. We tested the significance of this correlation 
using the Spearman's rank correlation method. We ob- 
tained a rank correlation coeflicient of 0.94 significant at 
the greater than 99.997% confidence level for the correla- 
tion between photon index and 0.2 2 keV luminosity and a 
coefficient of 0.71 significant at the 98.6% confidence level 
for the correlation between photon index and 2-10 keV 
luminosity. We note that the measured rest-frame 0.2- 
2 keV luminosity for high redshift quasars obtained from 
Chandra and XMM-Newton data is less accurate than the 
rest-frame 2-10 keV luminosity because it depends on ex- 
trapolation and suffers from the uncertainty in the absorb- 
ing columns towards the quasars. We adopt a significance 
at the 98.6% confidence level for this correlation, evalu- 
ated with the 2-10 keV luminosities. We also tested for a 
possible correlation between X-ray luminosity and redshift 
and between photon index and redshift and found none. 

The correlation between the photon index and X-ray lu- 
minosity was a surprising result. Therefore, we tested if 
this correlation is driven by certain data points. Exclud- 
ing one of our data points at a time and performing the 
Spearman's rank correlation between T and the 0.2-2 keV 
luminosity to the rest of the data set, wc find a strong cor- 
relation with a significance level greater than 99.98% each 
time. We performed a similar analysis with the 2-10 keV 
luminosity, and the correlation is significant at the greater 
than 95% confidence level each time. This indicates that 
the correlation is not driven by any particular data point. 
We further performed correlation tests by excluding data 
points with errors in F larger than 0.20. Two data points 
are excluded and the rest of the data points still show 
a correlation significant at the greater than 99.98% con- 
fidcnc;c level. The significance is at the 92% confidence 
level when 2-10 keV luminosities are used. We also tested 
if there is a bias from different flux levels in our sample 
such that quasars with low fluxes have a systematically 
flatter slope because of measurement errors. In principle, 
this should not be a problem for the sample in this pa- 
per because some low luminosity quasars have high S/N 
spectra as a result of the gravitational lensing effect. To 
test this, wc simiilatcd spectra with the same photon in- 
dex (r = 2) but with a wide range of S/N and then fitted 
them with the same models within XSPEC. We found that 
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for low S/N spectra, when the total number of photons 
is less than ^200, the measured F could possibly be flat- 
ter by ~ 0.1. This would afi^ect three of our data points: 
HE 2149-2745, LBQS 1009-0252, and Q 1208+101 with 
(one low-luminosity and two high- luminosity quasars) . Wc 
tested this effect by increasing the photon index by 0.1 
for these three quasars and performed the correlation test 
again and found the correlations between T and 0.2-2 
keV (2-10 keV) luminosities are significant at the 99.99% 
(98.7%) confidence levels, respectively. 

A T-Lx correlation has been previously searched for in 
other samples of radio-quiet quasars. Reeves et al. (1997) 
reported first a possible correlation between F and lumi- 
nosity with nine radio-quiet quasars observed with ASCA. 
However, this correlation was later not found in the study 
of a larger sample of 27 radio-quiet quasars observed with 
ASCA (including the original nine) by Reeves & Turner 
(2000). George et al. (2000) also searched for a cor- 
relation between F and luminosity in another sample of 
radio-quiet quasars observed with ASCA and did not find 
one. However, the quasars in our sample all have rela- 
tively high redshifts 1.7 < z < 4 compared to the red- 
shifts of quasars incorporated in previous studies of the 
T-Lx correlation. There is only one out of 26 quasars 
with z > 1.5 in the sample of George et al. (2000), and 
there are six out of 27 quasars with z > 1.5 in the sam- 
ple of Reeves & Turner (2000) and most of them have 
large errors on F. In addition, the high redshift quasars in 
the sample of George et al. (2000) and Reeves & Turner 
(2000) have high luminosities and span a small luminosity 
range of 10'*^-10'*^ erg s^^. The sample of lensed quasars 
in this paper have lower luminosities and span a larger lu- 
minosity range of 10^' 10^^ erg s~^. Recently, Page et al. 
(2003) presented another sample of quasars observed with 
XMM-Newton which contains several high redshift, low- 
Lx quasars. Wc plotted the T-Lx diagram for quasars of 
z > 1.5 from the sample of Page et al. (2003), Reeves & 
Turner (2000), George et al. (2000), and Vignah et al. 
(1999) in Figure 4. We excluded one data point from Page 
et al. (2003) with a large measurement error on F. We 
also excluded PHL 5200 from Reeves & Turner (2000) 
in our analysis because new XMM-Newton observations of 
the object showed that most of the X-rays originate from a 
nearby radio source (Brinkmann, Ferrero, & Gliozzi 2002). 
We corrected the X-ray luminosity of HE 1104—1805 re- 
ported in Reeves & Tmncr (2000) to accoimt for the grav- 
itational lensing magnification of about 12 based on our 
analysis presented in §4. We used rest-frame 2-10 keV 
luminosities in order to be consistent with the previous 
analysis. Most of the data from Page et al. (2003) are 
consistent with the T—Lx correlation found in our sam- 
ple of lensed quasars, especially in the low Lx range be- 
tween 10*^-^-10^^'^ erg s~^. On the high luminosity end 
(Lx > 3 X lO"*'"^ erg s'^) of high redshift quasars (z > 1.5), 
the F dependence on Lx seems to flatten out or even has 
an anti-correlation pattern. We performed the Spearman's 
rank correlation test to all of the {z > 1.5) quasars and did 
not find a correlation between F and 2-10 keV luminosity. 
We separated the quasars in two groups with luminosi- 
ties below and above 3 x lO"'^ erg s~^ and performed the 
rank correlation again to each group. The low luminosity 
quasars show a strong correlation between F and 2-10 keV 



luminosity at the 99.97% confidence level and the high lu- 
minosity quasars show an anti-correlation between F and 
2-10 keV luminosity at the 98.9% confidence level. 

We performed a Monte-Carlo simulation to test how 
much the F Lx correlation is affected by the uncertainty 
of the magnification factors of our lensed quasars. We used 
the combined data sets from the lensed sample and from 
the quasars of Page et al. (2003) with z > 1.5, which 
show a strong T-Lx correlation. We assumed a system- 
atic uncertainty on the magnification of a factor of two 
(see discussion in §4) and simulated 10,000 data sets with 
Lx perturbed from its measured value randomly, within 
the range of this error-bar. We fixed the data points from 
Page et al. (2003) at their original values. We performed 
the Spearman's correlation test between F and 2-10 keV 
luminosity on the simulated data sets and found that the 
T-Lx correlation in each simulated data set is at least 
significant at the 98.6% confidence level for 10,000 simula- 
tions and 9,888 of them have F Lx correlations significant 
at the greater than 99.5% confidence level. Therefore, it is 
unlikely that this correlation is driven by the uncertainties 
in the magnification factors. 

Possible errors in magnification factors, the small num- 
ber of quasars in the present sample combined with the 
medium S/N of several of the observations may have led 
to an unaccounted systematic effect. Additional observa- 
tions of 2 ~ 2 lensed quasars with better S/N and a more 
detailed modeling of all the lenses in our sample will help 
confirm this result. Although the T-Lx correlation found 
in this paper could be a result of selection effect, we discuss 
possible interpretations of this correlation below. 

Previous studies did not find any correlation between 
F and Lx in low-redshift quasars and Seyfert Is. Re- 
cently, several X-ray monitoring programs of Seyfert galax- 
ies have indicated a correlation between X-ray spectral in- 
dex and 2-10 keV X-ray flux (Chiang et al. 2000; Petrucci 
et al. 2000; Vaughan '& Edelson 2001). The spectra of 
these Seyfert galaxies steepen as the X-ray flux increases. 
Vaughan & Edelson (2001) suggested that variations of 
several properties of the X-ray emitting corona (e.g., elec- 
tron temperature, optical depth, size) with the X-ray flux 
may explain this F vs. X-ray-flux relation. Similarly, a F 
vs. X-ray-flux relation also exists in some Seyfert 2 and 
radio-loud objects (Georgantopoulos & Papadakis 2001; 
Zdziarski & Grandi 2001; Gliozzi, Sambruna, & Eracleous 
2003). These observations indicate that although the F 
Lx relation does not apply to a sample of low-redshift 
quasars and Seyferts, this relation applies to individual 
objects. In contrast to the low-redshift quasars, the F- 
Lx dependence found at high redshift, especially in high- 
redshift and \ow-Lx quasars indicates that these quasars 
may represent a distinct sub-population of quasars. It is 
possible that the properties of the X-ray emitting coronae 
of the high-redshift quasars are more homogeneous than 
those in the low-redshift quasars. 

6.2. Possible Interpretations of the Correlation Between 
the X-Ray Luminosity and Spectral Index 

We explore two possible interpretations of the correla- 
tion between the spectral index and the X-ray luminosity 
in the context of the model of Haardt & Maraschi (1993) 
and Haardt, Maraschi, & Ghisellini (1994). Our first hy- 
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pothesis is that there is a narrow range (of order a few) 
of accreting black hole masses at high redshift while the 
observed range of observed X-ray luminosities (spanning 
aproximately two orders of magnitude) is caused by a large 
range in the accretion rate, i.e., a wide range of Edding- 
ton ratios (the ratio of the accretion rate to the Eddington 
rate). Our second hypothesis is that the opposite is true, 
namely that the range of accreting black hole masses at 
high redshift is large (spanning approximately two orders 
of magnitude) , while the Eddington rate is fairly constant 
and close to unity (i.e., the accretion rate is always fairly 
close to the Eddington limit, regardless of the mass). 

In the model of Haardt & Maraschi (1993) and Haardt, 
Maraschi, & Ghisellini (1994) the inner accretion disk 
is sandwiched by a hot, tenuous, possibly patchy corona 
which emits the hard X-ray photons. The corona cools by 
inverse Compton scattering of soft photons {E < 100 eV) 
from the underlying disk, which results in the hard X-ray 
photons that we observe. Haardt, Maraschi, & Ghisellini 
(1997) explored how the model predictions change for a 
wide range of the system parameters, namely the optical 
depth and temperature of the corona and the temperature 
of the (assumed) black-body spectrum of seed photons. In 
particular, Haardt, Maraschi, & Ghisellini (1997) found 
that the 2-10 keV spectral index is nearly proportional to 
the optical depth in the hot corona and also decreases as 
the temperature of seed soft-photon spectrum increases. 
Moreover, in the case of a "compact" corona, i.e., one in 
which pair production is high enough that the Compton 
scattering optical depth is dominated by electron-positron 
pairs, the 2-10 keV spectral index is nearly proportional to 
the log of the 2-10 keV luminosity, just as we have found 
observationally. 

In our first hypothesis the observed range of luminosi- 
ties is set by the range of accretion rates, with the black 
hole mass being nearly constant. In this context the ob- 
served correlation might be explained if the optical depth 
in the hot corona is somehow related to the accretion rate; 
for example, the density of the corona could increase with 
accretion rate and so would the optical depth. There is a 
competing effect, however: as the accretion rate decreases 
so does the temperature of the underlying thin disk, which 
provides the seed photons, with the result that the spectral 
index of the emerging X-ray spectrum increases. 

In our second hypothesis, the accretion rate for all ob- 
jects is very close to the Eddington limit, while the ob- 
served range in luminosity is set by a range in black hole 
masses. In the context of this hypothesis, and under the 
assumption that the optical depth of the corona is dom- 
inated by electron-positron pairs, Haardt, Maraschi, & 
Ghisellini (1997) predict a correlation between the 2- 
10 keV spectral index and luminosity just like the one we 
observe. The optical depth in this case depends only on 
the "compactness" of the corona, which can be expressed 
as 4 ~ 10*{h/r^){Lx/LEdd) (Haardt & Maraschi 1993), 
where Lx/LEdd is the ratio of the X-ray luminosity to 
the Eddington luminosity (roughly proportional to the Ed- 
dington ratio) and h and r are the vertical and radial ex- 
tent of the corona in units of the Schwartzschild radius. 
Thus, in this scenario, the spectral index could increase 
if the compactness of the corona increases with black hole 
mass (See equation 16 of Haardt & Maraschi (1993)), with 



the consequence that the emerging 2-10 keV luminosity 
increases as well. Moreover, as the optical depth reaches 
unity, the spectral index reaches its maximum value and 
begins to decline for larger optical depths. This feature 
may explain the turnover that we observe in the spec- 
tral index vs luminosity plot of Figure 4. We therefore 
favor this interpretation over the previous one. This in- 
terpretation is also appealing because it is consistent with 
the predictions of semi-analytic models of Kauffmann & 
Hachnelt (2000) for the cosmological evolution of super- 
massive black holes and their fueling rates. In particular, 
these authors predict that at redshifts between 1.5 and 3, 
the fueling rate of supcrmassive black holes in quasars are 
within a factor of a few of the Eddington limit, while their 
masses span the range between lO'' and 10^ M©. 

6.3. The Optical-to-X-Ray Index, 

Vignali, Brandt, & Schneider (2003); Vignali et al. 
(2003a, b) found that the optical-to- X-ray index, aox, of 
radio-quiet quasars is mainly dependent on the rest-frame 
ultra-violet luminosity such that quasars with higher ultra- 
violet luminosities have steeper aox values (see Figure 5). 
They also do not detect a significant redshift dependence of 
aox in their analysis, however, Bechtold et al. (2003) found 
that aox depends primarily on redshift (aox steepens with 
increasing redshift) and weakly on luminosity. We com- 
pare the average aox, -1-70 (-1.66 without HE 2149-2745), 
of our sample with the average aox, -1.71, of Vignali et al. 
(2003b), and the aox values are consistent within the ob- 
served r.m.s. dispersion. The redshift range of our quasars 
is 1.7 < z < 4, lower than the sample of z > 4 quasars 
studied in Vignah et al. (2003b). However, the consis- 
tency between the aox values from the lensed sample in 
this paper and the sample of z > 4: quasars studied in 
Vignali et al. (2003b) does not rule out a possible small 
evolution of aox with redshift due to our small sample size. 
We plotted the aox values from our lensed sample against 
their redshift and rest-frame 2500A luminosity density in 
Figure 5(a) and (b), respectively. We performed the Spear- 
man's rank correlation test between aox and z and between 
aox and L ° , but found no significant correlation in 

2500 A' ° 

either case. Considering the large dispersion of the aox 
values from the best fit ao^-L ° relation of Vignali, 

Brandt, & Schneider (2003), it is not surprising that no 
significant correlation is found in our data set because of 
the small sample size. We over-plotted our aox data on 

the aox-i ° relation of Vignali. Brandt, & Schneider 

2500A " ' ' 

(2003) with a dispersion of 0.25 as indicated by the shaded 
region in Figure 5(b), and found that most of our data are 
consistent with the aox-i ° relation except one BAL 

° 2500A ^ 

QSO data point. 
We also searched for possible correlations between the 

aox values and other properties of the qiiasars such as F 
and Lx- The correlation results were strongly biased by 
the -2.1 aox value of HE 2149-2745. When we removed 
this data point, no significant correlation was detected. 

6.4. Short Time Scale Variability 

We used the upper limits of the excess variance ob- 
tained from light-curves of bin size 1000 seconds when 
comparing with results for local Seyferts and other high 
redshift quasars. This bin size corresponds to rest-frame 
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time scales ranging from 200 to 370 s for the quasars in 
our sample and is similar to the time scales used in other 
short time scale variability studies in Seyfert galaxies and 
quasars (e.g. Turner et al. 1999). We plotted the ex- 
cess variances of the high redshift lensed quasar sample 
together with the variances for Seyfert Is (Nandra et al. 
1997), NLSls (Leighly 1999), LLAGN (Ptak et al. 1998), 
and high redshift {z > 2) quasars (Manners, Almaini, & 
Lawrence 2002) in Figure 6. We fitted the Seyfert 1 points 
with a power law and fixed the power-law index to the 
same value reported in Nandra et al. (1997), allowing 
only the normalization parameter to be free. The power- 
law relation is shown as a solid line in Figure 6. All of 
the upper limits of the excess variances for the high red- 
shift lensed quasars of this paper are consistent with the 
Seyfert 1 a^^^-Lx correlation. When we compare the up- 
per limits with the NLSls data, three of the upper limits 
of our sample arc located below the cr^^g Lx relation for 
NLSls, which have about an order of magnitude larger ex- 
cess variance than the Seyfert Is. The consistency between 
the excess variance upper limits of our lensed sample and 
the Seyfert Is a'^^ns relation does not contradict the 
possible upturn in the value of the variability of z > 2 
quasars for quasars with luminosities of about 10''^erg s~^ 
as presented by Manners, Almaini, & Lawrence (2002). 
Most of the quasars in our sample have redshifts greater 
than 2, but the luminosity range of our sample is below 
the luminosity range of Manners, Almaini, & Lawrence 
(2002). 

The excess variance analysis that we discuss in this pa- 
per probes the shortest time-scale variability of radio-quiet 
quasars. There are two main interpretations that have 
been presented in previous studies to explain the excess 
variance and luminosity relation, (e.g., Nandra et al. 1997; 
Manners, Almaini, fc Lawrence 2002). First, the more 
luminous sources may be physically large and lack the 
short time-scale variability. Second, the more luminous 
sources may contain more independently flaring regions 
and these flares contribute less to the total flux for the 
more luminous sources compared with the low luminous 
sources. There are several rapid flares observed in our 
light-curve sample, in particular, such flares were detected 
in RX J0911.4-h0551, Q 2237+0305, and PG 1115+080 
(Chartas et al. 2001; Dai et al. 2003, 2004). We per- 
formed the Kolmogorov-Smirnov test to the light-curves of 
the individual images of the systems, RX .10911.4+0551, 
Q 2237+0305, and PG 1115+080, and found the light- 
curves are variable at the 99.8%, 97%, and 99.6% confi- 
dence levels, respectively. These observations support the 
second explanation where there are also flares in the high 
luminosity quasars, however, due to the quasar's high lu- 
minosity, flares contribute less to the total luminosity and 
thus luminous quasars have smaller excess variances. 

7. CONCLUSIONS 



We presented results from a mini-survey of relatively 
high redshift (1.7 < z < 4) gravitationally lensed radio- 
quiet quasars observed with the Chandra X-ray Observa- 
tory and XMM-Newton. We demonstrated how gravita- 
tional lensing can be used to study high-redshift quasars. 
Our main conclusions are as follows: 

1. Wo find a possible correlation between the spectral 
slope and X-ray luminosity of the gravitationally 
lensed quasar sample. The X-ray spectral slope 
steepens as the X-ray luminosity increases. The 
limited number of quasars in the present sample 
combined with the medium S/N of several of the 
observations and the systematic uncertainties from 
the lensing magnification modeling may have led to 
an unaccounted systematic effect. Additional ob- 
servations of 2; ~ 2 lensed quasars with better S/N 
and more detailed modeling of all the lenses in our 
sample will allow us to confirm this result. 

Such a correlation is not observed in nearby z <0.1 

quasars suggesting that quasars at redshifts near 
the peak of their number density may have differ- 
ent accretion properties than low redshift quasars. 
When we combined the data from other samples 
of radio-quiet quasars selecting quasars with red- 
shift {z > 1.5) together with the present lensed 
sample, the correlation is still significant, espe- 
cially in the low X-ray luminosity range between 
2^043-45.5 pj.g g-i_ ^Y^is correlation is confirmed 
by future studies, it could provide significant infor- 
mation on the emission mechanism and evolution of 
quasars. 

We suggest that this correlation can be understood 
if we hypothesize that the quasars in our sample are 
fueled at rates near their Eddington limit (consis- 
tent with recent models for quasar evolution) and 
that the optical depth of their X-ray emitting coro- 
nae increases with black hole mass. 

2. We did not find a strong correlation between the 
optical-to-X-ray spectral index, aox, on either red- 
shift or on UV luminosity in our small sample. How- 
ever, most of our data points are consistent with 
the aoK-L ° correlation of Vignali, Brandt, & 

2500j4 

Schneider (2003). 

3. Our estimated upper limits of X-ray variability of 
the relatively high redshift lensed quasars sample is 
consistent with the known correlation observed in 
Seyfert Is. 
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Table 1 



The Sample of Gravitational Lensed Radio-Quiet Quasars 











First Chandra Observation 


Second Chandra Observation 


XMM-Newton Observation 








Vjaiaciic ivH 


Date 


Exposure 


Date 


Exposure 


Date 


Exposure 


Quasars 


Rcdshift 


bal? 


(10 cm ) 




(ks) 




(ks) 




(ks) 


HE 0230-2130 


2.162 




2.3 


2000-10-14 


15 










HS 0818+1227 


3.115 




3.4 


2002-12-18 


20 










APM 08279+5255 


3.911 


BAL 


3.9 


2000-10-11 


10 


2002-02-24 


90 


2002-04-28 


100 


RX J0911.4+0551 


2.80 


mini-BAL 


3.7 


1999-11-03 


30 


2000-10-29 


10 


2001-11-02 


16 


LBQS 1009-0252 


2.74 




3.8 


2003-01-01 


10 










HE 1104-1805 


2.303 




4.6 


2000-06-10 


49 










PG 1115+080 


1.72 


mini-BAL 


3.5 


2000-06-02 


26 


2000-11-03 


10 


2001-11-25 


60 


Q 1208+101 


3.80 




1.7 


2003-03-02 


10 










H 1413+117 


2.55 


BAL 


1.8 


2000-04-19 


40 










HE 2149-2745 


2.033 


BAL 


2.3 


2000-11-18 


10 










Q 2237+0305 


1.695 




5.5 


2000-09-06 


30 


2001-12-08 


10 







^The Galactic A'^h is based on Dickey & Lockmaii (1990). 



Table 2 

Fits to the spectra of Gravitational Lensed Radio-quiet Quasars 



o 



Observed Properties^ 











Lx (0.2-2 keV) 


Lx (2-10 keV) 


fx (0.4-8 keV observed) 






Quasars 


Heasnirt 


Model 


r 


(erg s~*) 


(erg s~^) 


^erg cm s j 






HE 0230-2130 


2.162 


pow 


I Qo+0.08 


5.7 X lO""^ 


4.5 X 10*^ 


2.3 X 10-1^ 


1.09(39) 


0.33 


HS 0818-1-1227 


3.115 


pow 


1 cQ+0.10 


2.9 X 10"^ 


4.4 X 10*^ 


1.3 X 10-1=* 


0.90(22) 


0.60 


APM 082794-5255'' 


3.911 


absori* (pow+zgau+zgau) 


- ,70+0.04 

l-78Io.o4 


3.5 X 10"° 


3.7 X 10*° 


4.7 X IQ-i^ 


0.86(235) 


0.95 


APM 082794-5255" 


3.911 


absori* (pow+zgau+zgau) 


1 94+0.02 


5.3 X 10''° 


4.3 X 10*° 


5.0 X 10-1=* 


1.12(168) 


0.15 


RX J0911+0551 


2.8 


zwabs*(pow+zgau) 


, o,+0.09 


9.0 X lO** 


2.6 X 10*= 


1.1 X lO^i-'' 


0.99(26) 


0.47 


LBQS 1009-0252 


2.74 


pow 




4.5 X 10"*^ 


3.6 X 10*^ 


1.0 X lO"^-'' 


1.57(9) 


0.12 


HE 1104-1805 


2.303 


pow 


1 Sfi+006 


4.4 X lO"**^ 


4.0 X lO"-'' 


1.8 X 10"" 


0.97(82) 


0.55 


PG 1115-1-080 


1.72 


absori* (pow+zgau+zgau) 




8.1 X lO"*^ 


6.5 X 10*^ 


5.3 X 10"" 


1.15(200) 


0.07 


Q 1208+101 


3.8 


pow+zgau 


„ , + 0.2 
^•■'-0.2 


1.4 X 10*" 


6.8 X 10**^ 


7.1 X 10-" 


1.13(6) 


0.34 


H 1413+117 


2.55 


zwab*(pow+zgau) 


1 s + 0-3 

-■-■0-0.3 


4.0 X lO""^ 


4.5 X 10**^ 


1.0 X 10-1^ 


0.86(19) 


0.63 


HE 2149-2745* 


2.033 


zwabs*pow 


1 4+0-6 


1.1 X 10" 


2.6 X lO** 


2.0 X 10"" 






Q 2237+0305 


1.695 


zwabs* (pow+zgau) 


1 no+0.07 


6.5 X 10""^ 


5.9 X 10**^ 


4.6 X 10-" 


1.13(153) 


0.14 



Note. — The spectral fits were performed within the observed energy ranges 0.4—8 keV, and all derived errors are at the 68% confidence level. 

'^All models contain Galactic absorption model wabs, and all Chmidra spectra contain ACTSABS model lo correct the contamination from the ACTS filter. The 
pow, absori, zgau, and zwabs models represent powerlaw, ionized absorber, redshifted Gaussian emission or absorption feature, and redshiftcd absorption models, 

respectively. ^ 

^The luminosity bandpasses are given in the rest-frame of the quasar. The various absorption components are removed when calculating luminosites. The observed 
values have not been corrected for the lensing magnification. H 

^P(x^/^') is the probability of exceeding for v degrees of freedom. 

"^The fit results correspond to the Chandra spectrum of APM 08279+5255. ^ 
^The fit results correspond to the XMM-Newton spectrum of APM 08279+5255. 

^The C statistic is used when fitting this low S/N spectrum. Other Chandra and XMM-Newton spectra are binned with 15 and 100 events per bin, respectively. 
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Table 3 

aox Values of Gravitationally Lensed Quasars in the Sample 



Quasars 


Optical 
Information^ 


''2500A 


I 

2500A 


f d 
J2keV 




HE 0230-2130 


CASTLES 


0.37 + 0.07 


1.2 


11. + 1. 


-1.35 + 0.04 


HS 0818+1227 


castles 


1.9 ± 0.4 


8.4 


9. + 1. 


-1.65 + 0.05 


APM 08279+5255 


BOO" 


2.59 ± 0.07 


13. 


5.4 + 0.4 


-1.80 + 0.01 


RX J0911+0551 


CASTLES 


1.2 + 0.2 


4.9 


2.3 + 0.3 


-1.81 + 0.04 


LBQS 1009-0252 


castles 


3.7 + 0.7 


15. 


32. + 6. 


-1.56 + 0.05 


HE 1104-1805 


CASTLES 


3.0 ± 0.6 


10. 


11.0 + 0.6 


-1.70 + 0.03 


PG 1115+080 


C87' 


1.7 ± 0.1 


4.4 


11.9 + 0.4 


-1.59 + 0.01 


Q 1208+101 


CASTLES 


11. ± 2. 


54. 


60. + 20. 


-1.62 + 0.06 


H 1413+117 


CASTLES 


1.6 ± 0.3 


6.1 


5. + 3. 


-1.7 + 0.1 


HE 2149-2745 


CASTLES 


7. + 1. 


22. 


2. + 1. 


-2.1 + 0.1 


Q 2237+0305 


OGLE 


8.0 + 1. 


22. 


16. + 1. 


-1.81 + 0.02 



^Thc programs or papers where we obtained the optical magnitudes used in the a, 

calculation. 

^Rest-frame 2500 A flux density, in units of 10~^^ ergs cm~^ s~^ Hz~^. 

■^Rcst-framc 2500 A luminosity density, in units of 10"^^ ergs s""'" Hz""'". 
■^Rest-frame 2 keV flux density, in units of 10""^^ ergs cm~^ s~^ Hz~^. 
^Egami et al. (2000). 
^Christian, Crabtree, & Waddell (1987). 
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Table 4 

Magnification of Gravitationally Lensed Quasars in the Sample 



Quasars 


Literature 


Model 


Flux Magnification 


Lx (0.2 - 2 keV) 
(10** erg s-^) 


Lx (2 
(10** 


- 10 keV) 
erg s~^) 


HE 0230-2130 


602"* 




15 


3.8 


3.0 




HS 0818+1227 


B02 




10 


2.9 


4.4 




APM 08279+5255 




SIE 


142 


2.4 


2.6 




RX J0911+0551 




SIE 


17 


0.53 


1.5 




LBQS 1009-0252 




SIE 


3.5 


12.9 


10.4 




HE 1104-1805 




SIE 


12 


.3.6 


3.3 




PG 1115+080 


T02'' 


JF+CUSP+SIS" 


26 


3.1 


2.5 




Q 1208+101 


B02 




3.1 


45.8 


21.7 




H 1413+117 


C99' 




23 


1.7 


2.0 




HE 2149-2745 




SIE 


3.4 


0.33 


0.78 




Q 2237+0305 


898=' 




16 


4.0 


3.7 





"Schmidt, Webster, & Lewis (1998) 
'■Treu & Koopmans (2002) 

*^JF, CUSP, and SIS are JafFe, Cuspy, and Singular Isothermal Sphere mass models for the luminous mass of the 
galaxy, dark part of the galaxy, and nearby group of the system, respectively. 

''Barvainis & Ivison (2002) 

°Chae & Turnshek (1999) 
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Table 5 

Excess Variance of Gravitational Lensed Radio-quiet Quasars 



Quasars 


(10-2) 
Bin Size = 500 s 


(10-2) 
Bin Size = 1000 s 


(10-2) 
Bin Size = 2000 s 


HE 0230-2130 


-1.7 ± 0.7 


-0.5 + 0.5 


-0.4 + 0.3 


HS 0818+1227 


0.1 ± 2.2 


-0.6 + 2.0 


-0.3 + 0.8 


APM 08279+5255 


-0.5 ± 0.3 


-0.3 + 0.2 


-0.1 + 0.1 


RX J0911+0551 


-3.1 ± 2.3 


-0.7 + 2.3 


0.9 + 2.9 


LBQS 1009-0252 


-0.5 ± 2.6 


-3.1 + 0.8 


-1.4 + 0.9 


HE 1104-1805 


1.9 ± 1.1 


1.1 + 0.8 


1.0 + 0.8 


PG 1115+080 


-0.2 ± 0.5 


-0.2 + 0.3 


0.1 + 0.3 


Q 1208+101 


-6.1 ± 1.5 


-2.7 + 0.9 


-1.2 + 0.8 


H 1413+117 


-2.4 + 3.5 


-3.0 + 2.5 


0.0 + 2.3 


Q 2237+0305 


-0.6 + 0.3 


-0.3 + 0.1 


0.0 + 0.2 
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Table 6 



Excess Variance of Quasars Observed with XMM 



Quasars 


(10-5) 
Bin Size = 500 s 


2 

(10-5) 
Bin Size = 1000 s 


2 

(10-5) 
Bin Size = 2000 s 


2 

(10-5) 
Bin Size = 3000 s 


APM 08279+5255 


0.4 ± 0.2 


0.1 ± 0.2 


0.1 ± 0.2 


0.1 ± 0.1 


RX J0911+0551 


4.8 ± 1.5 


3.2 ± 1.5 


1.4 ± 1.0 


-0.1 ± 0.4 


PG 1115+080 


0.2 ± 0.2 


0.1 ± 0.2 


0.0 ± 0.1 


-0.11 ± 0.03 
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Fig. 1. — Spectra of gravitationally Icnscd radio-quiet quasars observed with Chandra, except for the spectra of APM 08279+5255 and 
PG 1115+080 that were observed with XMM-Newton. (Part I) 
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Fig. 1. (continued) 
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Fig. 2. — Light-curves of gravitationally lensed radio-quiet quasars observed with Chandra. The light-curves are binned with a bin size of 
1000 s in the observed-frame. 
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Fig. 3. — Plots of (a) X-ray Luminosity vs. redshift. The squares represent 0.2-2 keV luminosities and crosses represent 2-10 IceV 
luminosities. (b)X-ray photon index vs. redshift. (c) X-ray photon index vs. 0.2-2 keV luminosity, (d) X-ray photon index vs. 2-10 keV 
luminosity. 
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(erg s^) [2-10 keV Restframe] 



Fig. 4. — T-Lx (2—10 keV rest-frame) diagram for high rcdshift quasars (z > 1.5). The data shown as filled triangles are from Page et al. 
(2003). The data shown as filled circles arc from Reeves & Turner (2000). The data shown as crosses are from George et al. (2000). The 
data shown as open circles are from Vignali et al. (1999). The data shown as open squares are from our own sample of lensed quasars. 
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Fig. 5. — Plots of Qox vs. redshift and UV luminosity. The crosses represent BAL QSOs and squares represent non-BAL quasars for all 
of the plots in this figure, (a) Oox vs. redshift. (b) aox vs. rest-frame 2500A luminosity. The solid line represents the best fit ciox-L ° 
relation obtained from Vignali, Brandt, & Schneider (2003), and the shaded region represents a dispersion of 0.25 from the best fit line. 
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Fig. 6. — Excess variance in the 0.2-10 keV band vs. 2-10 keV luminosity. The data points representing the Seyfert Is shown as filled 
triangles were obtained from Nandra et al. (1997) and are fitted with a power law (solid line). The ASCA SIS light-curves of the Seyfert 
Is were binned in 128 s bins. The data points for the NLSls shown as open circles were obtained from Leighly (1999). The ASCA SIS 
light-curves of the NLSls were binned in 128 s bins. The data points for the LINERS and LLAGNs shown as squares were obtained from 
Ptalc et al. (1998). The data points shown as filled stars were obtained from z > 2 quasars of Manners, Almaini, & Lawrence (2002). The 
downward arrows represent upper limits for the excess variance of the gravitationally lensed quasars from the present sample. 



